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We show how the scattering phase shift, the s-wave scattering length and the p-wave scattering 
volume can be obtained from Riccati equations derived in variable phase theory. We find general 
expressions that provide upper and lower bounds for the scattering length and the scattering volume. 
We show how, in the framework of the variable phase method, Levinson's theorem yields the number 
of bound states supported by a potential. We report new results from a study of the heteronuclear 
alkali dimers NaRb and LiRb. We consider ab initio molecular potentials for the X 1 E + and a 3 E + 
states of both dimers and compare and discuss results obtained from experimentally based X 1 £ + 
and a 3 E + potentials of NaRb. We explore the mass dependence of the scattering data by considering 
all isotopomers and we calculate the numbers of bound states supported by the molecular potentials 
for each isotopomer. 
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I. INTRODUCTION 

The study of the properties of ultra-cold trapped en- 
sembles of atoms requires a description of very slow 
atomic collisions. The scattering length is a crucial quan- 
tity in the study of Bose-Einstein condensation. While 
there has been much interest in the behaviour of ultra- 
cold monatomic alkali gases, recent successes in the pro- 
duction of Bose- Fermi ensembles of various heteronuclear 
alkali dimers EJ 0, Q and dual species Bose-Einstein 
condensates |4j have created further interest in studies 
of binary mixtures. Scattering of heteronuclear atoms 
at ultra-low temperature is characterised by the s-wave 
scattering length, a s , the effective range, i? e ff, and the p- 
wave scattering volume, a p , the contributions from higher 
angular momenta being negligible at very low energy 
H H Q, H 0]. There is a need for the values of the 
scattering lengths and scattering volumes. 

In this paper we make use of the variable phase ap- 
proach to potential scattering 0] to calculate the scat- 
tering data a s and a p . This method is convenient and 
to some extent superior to the traditional calculation 
of phase shifts. It is advantageous compared to the 
usual calculations of the solutions of a second-order lin- 
ear differential equation, the radial Schrodinger equation, 
that are matched to combinations of asymptotic solu- 
tions at large atomic separations, R. In the variable 
phase method one solves a first-order nonlinear differen- 
tial equation, the phase equation, of the Riccati type, 
whose solution in the asymptotic region provides the 
phase shift directly. As we show below, the phase equa- 
tion is simple to manipulate in deriving properties of the 
scattering. 

We have shown how to use the variable phase method 
to compute the s-wave scattering length of a pair of col- 
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liding atoms |ll| . We derived corrections arising from 
the long range interaction, accurate to at least first or- 
der in the interaction strength, and we provided upper 
and lower bounds to the scattering length. Use of the 
corrections significantly accelerates the convergence, as 
R — ► co, to the desired scattering length, a, of a quantity 
a(R) that is the scattering length accumulated at separa- 
tion R. Here we extend our previous work to any angular 
momentum. Our analysis is consistent with the exten- 
sion of the technique introduced for s- waves by Marinescu 
[l2j| . While we suggest that the the Riccati equation be 
solved directly, our corrections can also be applied to 
scattering parameters obtained from explicit computa- 
tion of the wavefunction. 

The variable phase method was demonstrated for s- 
waves 0| for the model potential described in reference 
13] . Here we compute the s-wave scattering lengths and 
p-wave scattering volumes for the heteronuclear dimers 
NaRb and LiRb. We use ab initio molecular potentials 
for the two states X X E + and a 3 E + as well as the 
latest X 1 E + and a 3 E + NaRb potentials based on exper- 
imental data [TEl fill . A difficult problem in the accu- 
rate determination of the scattering length arises from 
the quality of the interatomic potential. If the potential 
is not accurately known then it is possible that, within 
its error bounds, it supports a zero-energy bound state 
which makes the computed scattering length extremely 
sensitive to any change in the potential. Our results, ob- 
tained with the above potentials, are compared to other 
recent calculations fil| . 

We calculate the number of bound states supported by 
those potentials using the simplest formulation of Levin- 
son's theorem |ToL 1 1 8| . linking the scattering phase shift, 
Sk,i, to the number of bound states, N^, via N^tt — 
5k,i |fc=o where k is the wavenumber of the relative motion. 
The scattering phase shift may be computed from the 
phase equation of variable phase theory with no mod[7r] 
ambiguity 

This paper is organised as follows. In section 2 we 
present our mathematical model, recall results from vari- 
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able phase theory 0, 0] , give detail of our derivations 
and show that our results may also be derived by the 
method of Marinescu 0]. I n section 3 we discuss the 
construction of the molecular potentials [Til ITU Il6| and 
our results which we tabulate for all relevant isotopomers. 



coupled to Eq. with I = 0, is not amenable to numer- 
ical solution because of the poles that it contains; the 
variable phase method is not suitable for numerical cal- 
culation of the effective range but it can be used to find 
long range corrections to it. 



II. MATHEMATICAL MODEL 
A. The phase equation 

The traditional method to compute the scattering 
phase shift, 5k where the angular momentum in terms 
of the rationalised Planck's constant h is ITi, needs the 
solution of the Schrodinger equation at large separation. 
In the variable phase approach ^(j the phase shift is ob- 
tained directly from the accumulated phase, 5k,i(R) [Tof . 
The accumulated s-wavc phase satisfies the equation 



dR 



Sk,o(R) = -k^ViR) sin 2 [kR + 5 k>0 (R)] , (1) 



where V(R) = 2^V(R)/h 2 ) V(R) being the interaction 
potential and [i the reduced mass. The phase shift is the 
limit of the accumulated phase as R — > oo; it suffers no 
mod[7r] ambiguity. 



B. Low-energy scattering 

At angular momentum lh the tangent of the phase 
5k i if) can be expressed in effective range theory by 



C. Corrections to the calculated scattering 
parameters 

The function ai(R) converges rather slowly with in- 
creasing separation, R. We need the value at infinite 
interatomic separation, a;(oo). In practice we must stop 
the numerical calculation of the solution of Eq. © at a 
finite value, R Cl of R. Computation time and accumu- 
lated truncation error can be large because R c must be 
chosen to be sufficiently large that the long range part 
of the interatomic potential has neglible effect at sepa- 
rations greater than R r . Adopting the same analysis as 
in our previous work we extend our calculations to 
any value of the angular momentum lh and show how the 
Riccati equation, Eq. ©, can be used to derive correc- 
tions to be applied to the calculated scattering parame- 
ters, that compensate for stopping the calculation at the 
finite distance R c . 

We need the contribution of the long range interaction 
over the range [R c , oo] . The presentation is simpler in 
terms of dimensionless quantities. With the substitution 
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ai{R) =ai(R)/B l R\ 
where Bi = /ctj and the change of variable 



(4) 



tan S k ,i(R) 



■„2l + l 



[ai(R) + k 2 h{R) + 0(fc 4 )] , (2) 



in which the coefficients ai(R) satisfy the differential 
equation of Riccati type 



dR 



ai{R) 



J+) R i+i 



H 



a^R-'aiiR) V(R) (3) 



Z = (R/R ( 



\2l + l 



Eq. becomes 



where 



^ ~a l ^(Z-~a l ) 2 F l {Z), 



(5) 



(6) 



where the coefficients a\ + ^ and a\ ' denote 1 / (2/ + I)!! 
and (21 — I)!! respectively where (n + 2)!! = (n + 2)nU 
with 0!! = 1 = 1!!. The accumulated phase, Sk,i(R), is 
the phase shift that would be determined for a poten- 
tial truncated at finite separation, R, and hence Eq. J5J, 
which applies strictly for scattering by a short range po- 
tential, is valid. 

We arc interested in the coefficients ai(R). For / = 
and 1 = 1, and R — > oo, these coefficients are the s- 
wave scattering length, a s , and the p-wave scattering vol- 
ume, a p , respectively. Similar equations can be found for 
bi(R). In the asymptotic region the coefficient bo(R) is 
related to the effective range, R s, via bo(oo) = a 2 R c g/2. 
Unfortunately the equation satisfied by bo(R), which is 



Fl{Z) - {21 + lfZ 2 ■ 



(7) 



Eq. © has the same form as Eq. © of Ref. 01 with 
Fi(Z) in place of V(R). Hence, from Ref. we find 
that for an attractive dispersion potential and sufficiently 

large R c , upper and effective lower bounds, d[ U ^ and a\ L ^ 
respectively, are given by 



~(U) 

a l 
and 



a-L, 



-(1 - a^f W l<c +2 (I - S I>C ) X l<c -2Y l>c , (8) 
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~(L) 



- (1 - a hc f W t , c + 2(1- a IlC ) X Lc - 2Y lt , 



where 



and 



l + Xi >c -{l-ai, c )W lt , 



Wi{Z) = J Fi{Z)dZ, 
Xi{Z) = J Wi{Z)dZ, 



Yi{Z) = J X^dZ, 



(9) 



(10) 



(11) 



(12) 



with the subscript c denoting quantities evaluated at R = 
R c or Z = 1. Furthermore aj jC is an upper bound. 

Expressions in Eqs. Q, I|1U|) . I|ll|) and l|12|l are incon- 
venient for computation and, changing the variable of 
integration back to R, we find 



Wi(Z) = 



1 



21 + 1 



R. 



2/+1 



R~ 2l V(R)dR, 



(13) 



Xi(Z) = (2Z + 1) R ( 



-2/-1 



iT ( W((-R)d-R, (14) 



and 



= (2/ + 1) fl" 1 / R 2l Xi(R)dR. 



(15) 



Eqs. ©, ©, 03 ) OH) an d 03 can be used to deter- 
mine the scaled upper and lower bounds and the scaling 
removed, if necessary, by Eq. Q . 

Equation @ can also be determined by extending the 
method of Marinescu (l2^ to higher angular momenta. 
The radial Schrodinger equation is, at zero energy, 



The regular and irregular solutions of Eq. (|16|l , to first 
order in the potential strength, are 



Ul (R) = u\ 0) (R) [1 + B l R 2 c l+1 X l (R)]~2vl 0) (R)B l R 2 c l+1 Y l (R), 



and 



v l (R) = v?\R) [l + BtRf^XiiR)], 



(17) 



(18) 



where u\ '(R) = a\ + > R l+1 and vf'(R) = a\ ' R~ l are 
the regular and irregular solutions of Eq. in the ab- 
sence of the potential. Let yi :C be the solution of Eq. i|ltj|) 
at R = R c . The scattering parameter, a/. c , is given by 



(0)/ / (o) 

yi, c u]J -yi tC u\j 

(0)/ / (0) ' 



(19) 



and the scattering length, a9~\ to at least first or- 
der in the potential strength is given by 



(1) _ Rc^c - Vl,c u hc 



yi,cV[. 



VLcVl,c 



(20) 



where the prime denotes the first derivative with respect 
to R. From Eqs. I|17|) - (|20|l we find that the scaled quan- 
tity a { p / BiRf +1 can be identified with a\ L) of Eq. ©. 



D. Inverse power potential 



For the attractive inverse power potential 



V(R) 



-C n R n , 



(21) 



~ m (R) - m (R) - V(R) yi (R) = 0. (16) 



dR 



J 



there are restrictions on n and I; i?( 2 ' +1 ) (" 2 ) must van- 
ish as R — > oo. We find the lower bound 



(a«) 


1 


2a; iC , 5 2 c 


21 + 1 


(n- 2) -(2/ + 1) n-2 (n-2) + 


(22 + 1) 


1 




1 a; jC 




2Z + 1 


n-2 (n-2) + (2Z + l)_ 
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-2\V(n- 



1S 

•2) 



where 
(2 f iC n /ti 

the upper bound, 
correction in Eq. I|22[l 
the numerator of Eq. 



the dimensionless quantity 
RZ 1 , and the correction that 



and 



gives 

, is the numerator part of the 
The correction in the form of 
122|) was derived in a different 



way by Szmytkowski |2l| (whose scattering length is 
S; -1 x ai(oo)). 

There are remaining corrections £ ( L J and to be 
made to et[ L ^ and dj [ . They can be derived from 
Eqs. (J20J) and (|21I) of Ref. 03; as noted in Ref. P3 there 
are terms of third order in the potential strength to be 
added to the expression in Eq i|22[l . The upper bound, 
Eq. (JHJ, is correct to first order in the potential strength 
while the lower bound, Eq. 10, includes some second or- 
der terms; both need second and higher order corrections. 
For the potential of Eq. I|21l) the ratio of the upper and 
lower correction is given to first order asymptotically as 
R — > co by 



from which we can find an improved approximation 



(23) 



a i - 



p-1 

E. P-wave scattering 
For I = 1 Eq. © is 



(24) 



d3 



R~ 2 V(R). (25) 



whose solution, ai(i?), at infinite separation is the p-wave 
scattering volume. 

The improved approximation is, from eq. 1|24JI . 



3ap } + (2n - 5)oi L) 



1 2n-2 
S-waves are discussed in Ref. 



(26) 



F. Numerical method 

Direct finite difference methods cannot be used to solve 
Eq. © because the solution contains poles that corre- 
spond to the bound states supported by the potential. 
The equation can be solved by the log-derivative method 
|22ll23Ll24| or by changing the variables. Substituting 



R = tan[0(i?)], 



(28) 



in Eq. J2J we obtain 



d6>(0) _ cos 2 (6>) 
d</> sin 2 ((f)) 



tan 3 (</>) 



i 2 



tan(6») 



V[tan($], (29) 



a(R) = tan[6(R)], 



(27) 



which we can solve over a range [0, <j) c ] by the Runge- 
Kutta method [25| where 4> c , = arctani? c , is close to 
tt/2. 



G. Advantages of the variable phase method 

We have seen that the scattering parameter a; can be 
evaluated by finding a; jC (= ai(R c )) by propagating the 
solution of the first order differential equation, Eq. to 
a suitable interatomic separation R c and then applying 
the correction of Eqs. © and The variable phase 
method yields a simple derivation of these corrections. 
The corrections can also be applied to the value of az iC 
obtained by any other suitable method. 

An advantage of the variable phase method over the 
traditional method is that one needs to solve a differen- 
tial equation only once. In the traditional method one 
must solve a differential equation several times at dif- 
ferent small wavenumbers, extract the phase shifts from 
the asymptotic solutions, and use the extrapolation of 
Eq. J2|. Some experimentation is needed when choosing 
suitable wavenumbers. 

An alternative to this method described above is to 
solve the zero energy Eq. and extract the parameter 
ai from its solution by Eq. I|19|) . as was done by Mari- 
nescu |l2j| for ao- Here and in the variable phase method 
one must solve one differential equation once. However, 
because of rapid variation of the local wave number over 
the potential well a variable step of integration should be 
used to economise on computation and minimise accu- 
mulated truncation error. Adaptive methods, that select 
the step automatically, are well established for first or- 
der equations such as Eq. I|29|l . derived from Eq. 10 • 
Eq. is a second order equation. Automatic step con- 
trol is not so well developed for second order equations as 
for first order one and in using numerical methods such 
as that of Numerov, the step and points at which it is 
to be changed must be inserted explicitly. Thus the vari- 
able phase method is also advantageous over the method 
based on the zero energy equation. 

It should be noted that the log-derivative method of 
solving the variable phase Eq. J2J 11] does not have an 
advantage over the Numerov method for solving the zero 
energy equation; the considerations concerning explicit 
insertion of the step apply to both and also to the tradi- 
tional method described above. 



III. NUMERICAL RESULTS FOR NaRb AND 
LiRb 

A. The molecular potentials 

We compute scattering parameters for LiRb and 
NaRb. We use various potentials and compare and dis- 
cuss results. Korek et al. |14| calculated ab initio molec- 
ular potentials for the X 1 ^" 1 " and a 3 £ + states of LiRb 
and NaRb. Docenko et al. constructed an X 1 !] 4 " 

NaRb potential from experimental data and from the 



3 E + state potential made by Zemke et 
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In 



the long range dispersion potentials we used the van der 
Waals coefficients computed by Derevianko et al. 1261 |27l 
with Cq from Ref. [26| and Cs and Cio from Ref. |27j . 

There is difficulty in matching an analytic long range 
potential to the short range part which is usually avail- 
able only in tabulated form. Any artificial discontinuity 
at the matching point can greatly affect the calculated 
scattering data. For LiRb the matching point that gives 
the smoothest fit is at a separation of 40 bohr. We found 
that within the precision of the data of Korek et al. 0] 
the singlet and triplet potentials do not differ beyond 30 
bohr and hence we neglected exchange beyond the match- 
ing point at 40 bohr, which is outside the typical LeRoy 
radius !2gj of 10-15 A of alkali dimers. For NaRb we 
used the short range potentials of Korek et al. 0] and 
matched at 30 bohr to the long range dispersion form 
evaluated with the dispersion coefficients of Derevianko 
et al. [US- 

The computed scattering parameters crucially depend 
on the potential. It is important to compare results ob- 
tained with various molecular potentials. The most re- 
cent short range LiRb potential is that of Korek et al. 
[l4| . However Docenko et al. 0| provided new mea- 
surements pertinent to the X 1 E + molecular state of the 
NaRb molecule and Zemke et al. 0] also made calcula- 
tions for the a 3 E + state of the NaRb molecule. We use 
calculations with these potentials in a comparison with 
the recent results of Weiss et al. 0] . 

We include exchange explicitly only for the a 3 £ + state 
of the NaRb molecule for which the point for matching 
to the long range potential is within the LeRoy radius. 
We represent exchange by the expression of Smirnov and 
Chibisov [2^ in the formulation given by Weiss et al 0] . 



B. Numerical results 



We tabulate the calculated scattering length data, a s 
and a p , and the numbers of bound states, iVb, supported 
by the molecular potentials. The values are those ob- 
tained from the analyses above including corrections for 
a p , those that we obtained from our previous analyses 
for a s 0] and those obtained from Levinson's theorem 
applied to the solution of Eq. for N^. We compare 
our results with those for NaRb of Weiss et al. ■ 
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FIG. 1: Accumulated and corrected scattering lengths for 
6 Li 85 Rb in the triplet state [a 3 E + ] , as functions of the inter- 
atomic distance R. 



1. Study of LiRb 

We used the following values for the masses of 
the isotopes, all in atomic mass, units obtained from 
the National Institute of Standards and Technol- 
ogy (NIST): Me Li = 6.0151223, Mv Li = 7.016004, 
M S5Rb = 84.9117893 and M S 7 Rb = 86.9091835. For 
the singlet state X 1 S + we see that the nature of the 
isotopomer influences the value of the scattering length 
and its sign; the scattering length is positive for the 
smaller mass and negative for the bigger mass. However 
a change in the Rb mass changes the value of the scat- 
tering length but not its sign. The order of magnitude 
of the various values of the scattering length does not 
change. The greater influence of mass changes of Li 
compared to Rb is caused by the greater sensitivity of 
the reduced mass of the dimer to changes in the mass of 
the lighter atom. The scattering volume is more stable 
with its sign remaining unchanged and its values having 
the same order of magnitude and is less sensitive to 
a change of the masses of the atoms. The number of 
bound states is sensitive to the Li mass. 

The a 3 £ + potential well is much shallower than the 
X 1 S + well. Therefore the value of the reduced mass of 
the LiRb molecule does not dramatically affect the scat- 
tering length (although it becomes slightly negative for 
the heaviest isotopomer). The scattering volume keeps 
its sign and there is a change of only one bound state 
when the LiRb dimer is composed of either of the Rb iso- 
topes and the heavier of the Li isotopes. However there 
is a qualitatively opposite behaviour in the changes of a p 
when the mass of Li is greater as compared with the re- 
sults obtained for the X 1 !]" 1 " state; the scattering volume 
is greater in magnitude for the smaller values of the mass 
of Li in the X : E + state and smaller for the greater value 
of the mass of Li in the a 3 S + state. 

The semiclassical analysis of Gribakin and Flambaum 
|13| can be invoked to explain the variation of scattering 
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FIG. 2: Accumulated and corrected scattering volumes for 
6 Li 85 Rb in the triplet state [a 3 E + ] , as functions of the inter- 
atomic distance R. 
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FIG. 3: Accumulated phase shifts (in units of n radians) as 
functions of the interatomic distance R and number of bound 
states for 6 Li 85 Rb in the singlet and triplet states. 



length with mass. The semiclassical scattering length 
depends only on the action integral which, for a given 
potential, is proportional to the square root of the re- 
duced mass. The semiclassical formulas for the scattering 
length and the number of bound states, combined with 
the dispersion parameters of the potential and the com- 
puted scattering length, allow us to find the action inte- 
gral appropriate to our quantum calculations. The action 
integral divided by the square root of the reduced mass 
thus calculated should be independent of isotopomer. 
This quantity obtained from our quantal results varies 
by less than 2% over the range of LiRb isotopomers con- 
sidered, which is consistent with semiclassical theory. It 
has the values 1.6 and 0.49 atomic units for the X 1 E + 
and the a 3 E + potentials respectively. 

Figures 1 and 2 illustrate the convergence of the ap- 
proximations to the scattering parameters. Figure 3 
shows the accumulated phase shift, demonstrating how 
it is related to the number of bound states. 

New ab initio or spectroscopically determined LiRb 
molecular potentials for the X 1 S + and a 3 S + states will 
be extremely useful to allow comparison of the results 
displayed in Table [I] and hence to confirm the physical 
implications of the signs of the scattering parameters, 
and the numbers of bound states. This is important for 
NaRb; when we compare values of a s , a p and com- 
puted from three different potentials for each molecular 
state we find discrepancies. 



2. Study of NaRb 

We took the value for the mass of 23 Na given by NIST 
as 22.98976967 atomic mass units. There are two iso- 
topomers. 

Our results, displayed in Table ITT1 computed with the 
X 1 S + potential of Docenko et al. Ha are in excellent 
agreement with those of Weiss et al. 17] for 23 Na 85 Rb. 
Weiss et al. 01 claimed that for NaRb the ab initio 



TABLE I: The scattering length (in bohr), scattering volume 
(in bohr 3 ) and number of bound states for LiRb. 



LiRb 


A 1 E+ 


a B 


a p 


iVb 


6 Li 85 Rb 


-40 


-97720 


50 


6 Li 87 Rb 


-64 


-105899 


50 


7 Li 85 Rb 


26 


-44698 


54 


7 Li 87 Rb 


18 


-59648 


54 


LiRb 


a 3 E+" 


a s 


a p 


iVb 


6 Li 85 Rb 


26 


-40309 


15 


6 Li 87 Rb 


24 


-43583 


15 


7 Li 85 Rb 


3 


-78665 


16 


7 Li 87 Rb 


-0.25 


-82200 


16 



potentials of Korek et al. [14J may not be appropriate 
for the calculation of scattering lengths as they do not 
agree with the potentials determined from spectroscopic 
measurments. 

The values of a s and a p depend again on the mass of 
the molecule. However the number of bound states is not 
influenced by the choice of NaRb isotopomer because the 
Rb isotopes differ little in mass. The sign of the scatter- 
ing length is not influenced by the choice of NaRb iso- 
topomer for the singlet potential but, for the triplet state, 
there is a qualitative difference between the results com- 
puted from the potential of Korek et al. Q (negative 
scattering length) and those obtained from the potentials 
of Weiss et al. an d of Zemke et al. (positive scat- 
tering lengths). The numbers of bound states computed 
from potentials in Refs. 0, 0] and 0] agree whereas 
we find fewer bound states supported by both triplet and 
singlet potentials of Ref. [T^J . This explains the change 
of sign of the scattering length. 

We compare the scattering volumes computed from the 
potentials of Korek et al. [14J with those computed from 
the potentials of Zemke et al. [l5|[. A striking differ- 
ence is the opposite signs for both singlet potentials. In 
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TABLE II: The scattering length (in bohr), and number of 
bound states for NaRb. 



NaRb 




a s a 


a s b 


a s c 


a s d 


N h a 


N h b 


N b c 


N h d 


23 Na 85 Rb 


174 


62 


167 


178 


83 


76 


83 


82 


23 Na 87 Rb 


84 


29 


55 


87 


83 


76 


83 


82 


NaRb 


a 3 E+" 


a s a 


a s » 


a s c 


a B d 


N h a 


N h b 


N h c 


N b d 


23 Na 85 Rb 




- 51 


59 


105 




18 




22 


23 Na 87 Rb 




-102 


51 


91 




18 




22 



"Calculated from potential in Ref. LUl. 
'Calculated from potential in Ref. Il4l . 
c Results from Ref. Il7l or other ref eren ce therein. 
^Calculated from potential in Ref. Il5l . 



TABLE III: 


The scattering 


volume (in 


bohr 3 ) 


of NaRb. 


NaRb 


A X E+ 


a p a 


a p b 


a p ° 


a p d 


23 Na 85 Rb 


-681980 


52815 




-663097 


23 Na 87 Rb 


373425 


-114509 




459310 


NaRb 


a 3 E+" 


a p a 


a p b 


a p ° 


a p d 


23 Na 85 Rb 




-235850 




3085301 


23 Na 87 Rb 




-266190 




615416 



"Calculated from potential in Ref. 16]. 
'Calculated from potential in Ref. Il4l . 
c Results from Ref. Il7l or other ref eren ce therein. 
^Calculated from potential in Ref. 1 1 . 



a semiclassical analysis similar to that described above 



for LiRb we found that the action integral divided by 
the square root of the reduced mass obtained from our 
quantal results varies by less than |% for the NaRb iso- 
topomers considered. In atomic units it has the values 
1.4 and 0.32 for the X 1 S+ and the a 3 S+ potentials of 
Korek et al. 0] respectively, the values 1.5 and 0.39 for 
the X 1 E + and the a 3 E + potentials of Zemke et al. [TsT ] 
respectively and the value 1.5 for the X 1 E + potential of 
Docenko et al. [Tr3 |. 

IV. CONCLUSION 

We have shown that the variable phase method is ap- 
propriate for the study of very slow atomic collisions as 
it allows direct calculation of scattering parameters and 
the number of bound states. The method is convenient. 
It let us derive corrections that form upper and lower 
bounds for a s and a p . The corrections are simple and we 
suggest that use of the first-order correction with a mod- 
erate value of R c is sufficient to evaluate a s and a p . We 
demonstrated the simplicity and efficiency of the phase 
method with our study of Li-Rb and Na-Rb collisions, 
and we reported new scattering data for these systems. 
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